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Abstract: 

In this paper we discuss the problem of small frequency approximation of the causal dissipative 
pressure wave model proposed in Kowar et al. (2010). We show that for appropriate situations 
the Green function G c of the causal wave model can be approximated by a noncausal Green 
function G P M that has frequencies only in the small frequency range [—M,M] (M < 1/tq, tq 
relaxation time) and obeys a power law. For such cases, the noncausal wave G^ contains partial 



waves propagating arbitrarily fast but the sum of the noncausal waves is small in the L — sense. 
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1. INTRODUCTION: FREQUENCY FRAMEWORK 

In the first section we summarize the general framework 
of frequency driven wave dissipation. 

From the mathematical point of view, in case of a homo- 
geneous and isotropic medium, dissipative pressure waves 
can be modeled by (cf. e.g. Kowar et al. (2010)) 



p(x,t) 



with 



(G* X)t /)(x,t) 
dp 



x e 



t e 



p(-,t) = and ^(.,t)=0 for t<0, 



(1) 



where G denotes a distribution (Green function), * x . t 
denotes the space-time convolution and / denotes a forcing 
term (source term) which models wave generation. We 
note that / is the same forcing term as in the absence of 
dissipation as long as (1) holds. (This is not true if p(x, t) 
is calculated for t > to (to > 0) from initial data at to-) 

1.1 The Green function 

The Green function can be modeled by 

e -/3,(|x|,r-1 



G(x,w) 



47T Ixl 



g o I 



x e 



(2) 



where 9?(/3*(|x|,<j)) > is even in uj and 3(/3*(|x|, w)) 
is odd in u. Here G(x, ui) denotes the Fourier transform 
of G(x, t) with respect to time t. The last two conditions 
guarantee that G(x, t) is real- valued. We focus on the cases 

/3*(x,o;) = a*(u>) l x l ■ (3) 

We call a* = a*(aj) and a = 5R(a*(w)) the attenuation- 
dispersion law and the attenuation law, respectively. 

1.2 The wave equation 

The above dissipative wave model satisfies the following 
integro- differential equation 



D 



l_dp 
c dt 



p(x,t) = -/(x,t), (4) 



where D* denotes the time convolution operator 
F{D*(g)}{u) = -^=a*(uj)T{g}(uj) 



/2tt 

defined for appropriate functions g = g(t). Here J~{g} 
denotes the Fourier transform of g with respect to time. 

More details about attenuation-dispersion laws and re- 
spective wave equations can be found in e.g. Nachman 
et al. (1990), Szabo (1995), Waters et al. (2000), Chen 
et al. (2004), Patch et al. (2006), Kelly et al. (2008) 
and Kowar et al. (2012). 

2. ATTENUATION-DISPERSION LAWS AND 
CAUSALITY 

2.1 What do we mean by causality? 

Let cp < oo denote the (constant) speed of the wave 
front of the wave G. Causality requires that the wave front 
initiated at the origin at time t = arrives at position 
x 7^ not before the time period T(\x\) — j^f is over. 
Mathematically this is equivalent to the following causality 
condition: for each ci > cf we have 



G ( x, i + - — - 

ci 







if 



t < and |x| ^ . (5) 



Various dissipative wave models are analysed with respect 
to causality in Kowar et al. (2012). In particular, it is 
shown that (3) implies cp = const, if cp < oo. See 
also Kelly et al. (2008). 

2.2 Attenuation- dispersion laws (7 € (1,2],) 

In this paper we consider the attenuation-dispersion laws 



ai(— iuj) 



c Vl + (-itqlj) 



7-1 



and 



COS(77T/2) 



(6) 



(7) 



for uj e M with 7 G (1,2], Co € (0, 00) and a\, To, ai > 0. 
a pl — a* (uj) is called the frequency power law. 

Let 

a r 7 " 1 a 

a\ = — —2 — I cos(7 7r/2)| and a 2 = — , (8) 
2 c c 

then it follows from 

(-iw) 7 = M 7 |cos (7 |) - i sin (7 |) sgn(w)} 

that 

9f(a°(a;)) w — ai tan (7 ^ sgn(w) |w| 7 — a 2 w 

for |t ix)| 7_1 << 1. We see that a* l (uj) is a good ap- 
proximation of a%(uj) for sufficiently small frequencies. For 
the case To = 10~ 6 ^s (liquid), the small frequency range 
condition |row| 7_1 < 0.1 is visualized in Table 1. 

Table 1. Visualization of small frequency range 
condition |r o;| 7_1 < 0.1 for t = 10~ 6 ^s. 

7 range bound M 

~Tl H < 10~ 4 MHz 1O- 4 MHz 

1.5 |w| < 10+ 4 MHz 10+ 4 MHz 

2.0 M < 10+ 5 MHz 10+ 5 MHz 



2.3 Two results on causality 

In Kowar et al. (2012), it is proven that the speed cp of 
the wave front of G(x, t) with a% defined as in (6) satisfies 
cf < Co and in Kowar (2010) it is proven that cp = c$. 
Moreover, it is shown that G(x, t) with a*' defined as in (7) 
with a 2 = has not a bounded wave front speed. The same 
result is true for a 2 7^ 0. 

Loosely speaking, causality restricts the growth of a(uj) 
and consequently uj i-> G(x, uj) must not decrease too fast 
for fixed x € R 3 \{0}. If w h> G(x, w) decreases too fast, 
then it behaves like the truncated Green function 

G M (x,w) := G(x,w)x[_m,m]M, (9) 



for which causality condition (5) cannot hold. Here 

, , _ / 1 if uj e [-M, M] 
X[-M,M](w) :- otherwise 



(10) 



denotes the characteristic function of uj on [— M, M}. 

3. SMALL FREQUENCY APPROXIMATION 

We now derive two theorems that permit to estimate the 
quality of small frequency approximations. Let 

• G c (x, t) denote the Green function with a % defined as 
in (6), 



• G pl (x, t) denote the Green function with a% 1 defined 
as in (7), 

• Gm(x, uj) be defined as in (9) and 
. G M (x,t) := J-^GmXx,*). 

Theorem 1. Let M > and A , Ai, A3 > be such that 

a(M) +A \uj- M\ < a c {uj) < A x uj 2 + A 2 \uj\ (11) 
holds for \uj\ > M. For |x| ^ it follows that 

Li e v ' 



|G c (x,-)-G^(x,.)|| L2 < ,/2A 



|G c (x,-) 



|L 2 



X 



Proof, uj 1 ^ G c (x, w) and cj i->- G m (x, uj) are square 
integrable and thus their inverse Fourier transforms t i->- 
G c (x, t) and i i->- G c M (x,t) are square integrable, too. 
Because of the Plancherel-Parseval equality, it follows 



|G c (x,.)-G| f (x,.)|U 2 _ ||G c (x,.)-G c m (x,.)||l 2 



l|G c (x,.)|| L2 
From (11) and 



|G c (x,.)||l 2 



.(12) 



/ 



exp{— w 2 /2} duj = V2tt , 



it follows 

||G c (x,.)-G^(x,.)||i 2 < 2e -2«(M)|x| 



-2 A (w-M) |x| 



M 



(47r|x|)s 



■ dw 



< 



and 



|G c (x,.)||| 2 > 



e -2a(M) |x| 
(4 7T |x|) 2 A |x| 

g -2(A 2 \u\+A x u 2 ) |x| 

(4^N) 2 



■ dw 



(4 7r|x|)2y2A7R 
The theorem follows from (12) and the last two results. 

Remark. Because uj i-> G c m (x, uj) vanishes on a non- 
empty interval, causality condition (5) is not satisfied for 
G%j. However, Theorem 1 shows that G^(x, t) is a good 
L — approximation of G c (x, t) if M is sufficiently large. 

Theorem 2. Let B x := $t(a p J - a%), B 2 := Q?(a? - a%), 

C:= 1 — 2 e -^iC-> l-l C0S (B 2 (uj) |x|) + e~ 2 M 
and M 5 e (0, M] be such that 

||G^(x, W )||| 2 = (l-5)||G c (x,a;)||| 2 

for some < 5 « 1. For |x| ^ it follows that 



(13) 



' )IIl2 ^ = *) A(ixi) + ^ 2 (ixi) 

l|G c (x, -)\\ L 2 

with Dj-(|x|) := max ue ; 3 G(|x|,cj) 2 for h := [-M 5 ,M 5 ] 
and/ 2 :=MUi. 



Proof. As in the proof of Theorem 1, it follows that 



G c M (x, •) - <(x, .)\\ L , ||G c M (x, •) - G^(x, .)\\ L , 



/' 



2 d W 

|G^(x,.)| 2 dc 



|G c M (x,-)IU 2 ||G^(x,-)|U2 
The theorem follows from the last result, (13) and 

I^m( X , — <^mC 

< A(W) 
+£>2(|x|) | |G c M (x,-)| 2 du;. 

|«|>M { 

Cf. Fig. 3 from the example in Section 4 for which 8 s=s 6 • 
10~ 4 , M 5 ^ 10 MHz, J>i(|x|) » 5.6- 1 0~ 13 and £> 2 (|x|) w 
1.04, i.e. ^(1 -<*) £>i(|x|) +<S£> 2 (|x|) - 0.025. 

4. NUMERICAL EXAMPLE 

We now present a numerical example. We choose the 
following parameter values 



7 w 1.66, 



ai 



co 
2 c ai 



0.15 — , 



cos( 77 r/2)| 



t w 10->s, 
138.08 Mifz 



(where we have used (8) and ai w 0.04344^) and 

a 2 = — w 920.55 — 

Co cm 

to model a medium similar to castor oil. (For the relaxation 
time To see Kinslcr ct al. (2000) and for a\ and 7 see Szabo 
(1995)). 

The differences of the attenuation laws and the phase 
speeds of the attenuation-dispersion laws a.% and a*' are 
visualized in Fig. 1. From the figure we see that a*' 
is a good approximation of the attenuation-dispersion 
law a% for the (practically measurable) frequency range 
[0, 60] MHz. In the following we choose as small frequency 
bound: 

M = 100 MHz « — . 

TO 

Small frequency approximation of causal model. For the 
estimation in Theorem 1 we choose the following constants 

da pl 

A a = 0.7 — — (M), Ai = ai and A 2 = a 2 



and obtain 

||G c (x,.)-G^(x,.)||l 2 



< 0.0828- 



-4.877-10 b -|x| 



l|G c (x,.)|| L2 

Therefore the frequency range [— 10 2 , 10 2 ] is reasonable for 
a small frequency approximation of G c for distances larger 
than 10 -6 cm. 

Small frequency approximation via power law. In Fig. 3 
we have visualized the functions Mo h-> 9{Mq) := 
||G^ o (x, -)\\l 2 and w h-> G(|x|,w) for |x| = 1cm. Ac- 
cording to this figure and Theorem 2, G^ is a good 
L 2 — approximation of G c for the case |x| = 1cm. In 
Table 2 we have listed the relative error 

||G c M (x,.)-G^(x,.)|U 2 



£m(|x|) := 



for different distances L. We see that G^(x, t) and 

Table 2. Visualization of error €m for different 
distances L (in cm). 



10- 



10" 



10" 



10 



|G^(x,.)||l 2 



e M 7.62 ■ 10~ 8 7.35 -10~ 5 4.46 ■ 10~ 4 7.13 ■ 10~ 5 

G^f(x, t) are good approximations of G c (x, •) as long as 
M > 100 MHz and |x| > 10 -6 cm . 

In particular, this means that the wave can be predicted 
for 

10~ 6 cm R 
t > = 6.67- l(rV*. 

co 

Remark on causality. Because the wave G c has the 
finite wave front speed cf — co and G^(x, •) is a good 
L 2 — approximation of G c (x, •), it follows that G^ has 
in "some sense" the same wave front speed. In "some 
sense" means that G^ contains partial waves propagating 
arbitrarily fast but their sum is small in the L 2 — sense. 

Remark on large frequencies. We note that the power 
attenuation law a pl (uj) increases much faster than a c (uj) 
for large frequencies and that the respective phase speed 
c pl (uj) has a singularity at uj 1 w 7.950959 • 10 6 MHz and 
at u> 2 = —u>i (cf. Fig 2), respectively. In contrast to the 
latter fact, the phase speed c c (uj) of the causal model has 
no singularity. 

5. CONCLUSIONS 

We showed that - for appropriate situations - the causal 
wave G c defined by (2), (3) and (6) can be approximated 
by the noncausal wave G^ defined by (9) (2), (3), (7) 
and (8). In words, the wave G pl contains partial waves 
propagating arbitrarily fast but the sum of the noncausal 
waves is small in the L 2 — sense. This result fits in with the 
results in Kelly et al. (2008). 
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Fig. 1. Numerical comparison of the attenuation laws and 
phase speeds. The dashed lines correspond to model 
a*(w) and the solid lines correspond to model a* l (uj). 
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Fig. 2. Numerical comparison of the attenuation laws 
and phase speeds for large frequencies The dashed 
lines correspond to model a^(us) and the solid lines 
correspond to model a* l (oj). 
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Fig. 3. Visualization of functions Mq g(Mo) 
||G^f o (x, -)\\l 2 and u ^ C(|x|,oj) for |x| = 1cm. 



